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Abstract 



Brillouin zones were introduced by Brillouin |Br| in the thirties to describe quantum 
mechanical properties of crystals, that is, in a lattice in M". They play an important role in 
solid-state physics. It was shown by Bieberbach ^| that Brillouin zones tile the underlying 
space and that each zone has the same area. We generalize the notion of Brillouin Zones 
to apply to an arbitrary discrete set in a proper metric space, and show that analogs of 
Bieberbach's results hold in this context. 

We then use these ideas to discuss focusing of geodesies in orbifolds of constant curvature. 
In the particular case of the Riemann surfaces H^/r(fc) (fc = 2, 3, or 5), we explicitly count 
the number of geodesies of length t that connect the point i to itself. 



1 Introduction 

In solid-state physics, the notion of Brihouin zones is used to describe the behavior of an electron 
in a perfect crystal. In a crystal, the atoms are often arranged in a lattice; for example, in NaCl, 
the sodium and chlorine atoms are arranged along the points of the simple cubic lattice Z^. 
If we pick a specific atom and call it the origin, its first Brillouin zone consists of the points 
in which are closer to the origin than to any other element of the lattice. This same zone 
can be constructed as follows: for each element a in the lattice, let Loa be the perpendicular 
bisecting plane of the line between and a (this plane is called a Bragg plane). The volume 
about the origin enclosed by these intersecting planes is the first Brillouin zone, 6i(0). This 
construction also allows us to define the higher Brillouin zones as well: a point x is in 6„ if the 
line connecting it to the origin crosses exactly n — 1 planes Loa, counted with multiplicity. 

This notion was introduced by Brillouin in the 1930s ( |Bi|j ) , and plays an important 

also 



role in solid-state theory (see, for example, [ AM | or [|Io2(| , also |ri[). The construction which 
gives rise to Brillouin zones is not limited to consideration of crystals, however. For example, in 
computational geometry, the notion of the Voronoi cell corresponds exactly to the first Brillouin 
zone described above (see ]PS| ). We shall also see below how, after suitable generalization, this 
construction coincides with the Dirichlet domain of Riemannian geometry, and in many cases, 
to the focal decomposition introduced in Pel| (see also |Pe3| ] ) . 
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With some slight hypotheses (see section |2|), we generahze the construction of Brillouin 
zones to any discrete set 5 in a path-connected, proper metric space X. We generalize the 
Bragg planes above as mediatrices, defined here. 

Definition 1.1 For a and b distinct points in S, define the mediatrix (also called the equidistant 
set or bisector) Lab of a and b as: 

Lab = {x £ X \ d{x, a) = d{x, b)} . 

Now choose a preferred point xq in S, and consider the collection of mediatrices {Lx^^s} g^s- 
These partition X into Brillouin zones as above: roughly, the n-th Brillouin zone Bn{xQ) con- 
sists of those points in X which are accessible from xq by crossing exactly n — 1 mediatrices. 



(There is some difficulty accounting for multiple crossings — see definition 2^ for a precise 
statement.) 

One basic property of the zones Bn is that they tile the space X: 

\^ Bn{xi) = X and Bn{xQ) H Bn{xi) is small . 

Here, with some extra hypotheses, "small" means of measure zero. Furthermore, again with 
some extra hypothesis, each zone Bn has the same area. (This property was "obvious" to 
Brillouin). Both results were proved by Bieberbach in in the case of a lattice in M^. 
Indeed, he proves (as we do) that each zone forms a fundamental set for the group action of 
the lattice. His arguments rely heavily on planar Euclidean geometry, although he remarks 
that his considerations work equally well in M*^ and can be extended to "groups of motions in 
non-Euclidean spaces". In |Jol], Jones proves these results for lattices in W^, as well as giving 



asymptotics for both the distance from Bn to the basepoint, and for the number of connected 
components of the interior of Bn- In section ^ we show that the tiling result holds for arbitrary 
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discrete sets in a metric space. If the discrete set is generated by a group of isometries, we show 



that each i?„ forms a fundamental set, and consequently all have the same area (see Prop. 2.11 ) 



We now discuss the relationship of Brillouin zones and focal decomposition of Riemannian 
manifolds. 

If xi{t) and X2{t) are two solutions of a second order differential equation with xi(0) = 
X2(0) and there is some T / so that xi(T) = X2{T), then the trajectories xi and X2 are 
said to focus at time T. One can ask how the number of trajectories which focus varies with 
the endpoint x{T) — this gives rise to the concept of a focal decomposition (originally called a 



sigma decomposition). This concept was introduced in [Pel] and has important applications in 



physics, for example when computing the semiclassical quantization using the Feynman path 
integral method (see [Pe3| ) . There is also a connection with the arithmetic of positive definite 
quadratic forms (see |Pe2| , [ [KP[| , and Pe3| ) . Brillouin zones have a similar connection with 



arithmetic, as can be seen in section ^. 

More specifically, consider the two-point boundary problem 

X = f{t,X,x), x{tQ) = Xq, x{tl) = Xi, XjtjXjXGW. 

Associated with this equation, there is a partition of into sets S^, where {xq, xi,tQ,ti) is in 
Sfc if there are exactly k solutions which connect (xo,to) to (xi,ti). This partition is the focal 



decomposition with respect to the boundary value problem. In |PT|, several explicit examples 



are worked out, in particular the fundamental example of the pendulum x = — sinx. Also, 



using results of Hironaka (|Hi]) and Hardt ( [ [Ha[| ) , the possibility of a general, analytic theory 
was pointed out. In particular, under very general hypotheses, the focal decomposition yields 
an analytic Whitney stratification. 

Later, in |lKP| , the idea of focal decomposition was approached in the context of geodesies 
of a Riemannian manifold M (in addition to a reformulation of the main theorem of [|PT|| ). Here, 
one takes a basepoint xq in the manifold M: two geodesies 71 and 72 focus at some point y £ M 
if 71(7") = y = ^2{T)- This gives rise to a decomposition of the tangent space of M at x into 
regions where the same number of geodesies focus. 

In order to study focusing of geodesies on an orbifold (M, g) with metric g via Brillouin 
zones, we do the following. Choose a base-point pq in M and construct the universal cover X, 
lifting pq to a point xq in X. Let 7 be a smooth curve in M with initial point po and endpoint 
p. Lift 7 to 7 in X with initial point xq. Its endpoint will be some x S 7r~^(p). The metric g 
on M is lifted to a metric ^ on X by setting g = 7r*g. Under the above conditions, the group 
G of deck transformations is discontinuous and so vr^^(po) C X is a discrete set. One can ask 
how many geodesies of length t there are which start at pQ end in p, or translated to (X, 7), 
this becomes: How many mediatrices L^q^s intersect at x, as s ranges over tt~^{po)? 

Notice that if the universal cover of M coincides with the tangent space TM^, the focal 
decomposition of [ [KP| ] and that given by Brillouin zones will be the same. If the universal 
cover and the tangent space are homeomorphic (as is the case for a manifold of constant nega- 
tive curvature), the two decompositions are not identical, but there is a clear correspondence. 
However, if the universal cover of the manifold is not homeomorphic to the tangent space at 
the base point, the focal decomposition and that given by constructing Brillouin zones in the 
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universal cover are completely different. For example, let M be S", and let x be any point in it. 
The focal decomposition with respect to x gives a collection of nested n — 1-spheres centered at 
X] on each of these infinitely many geodesies focus (each sphere is mapped by the exponential 
to either x or its antipodal point). Between the spheres are bands in which no focusing occurs. 
(See ||Pe3| ] ) . However, using the construction outlined in the previous paragraph gives a very 
different result. Since is simply connected, it is its own universal cover. There is only one 
point in our discrete set, and so the entire sphere §" is in the first zone Bi. 

The organization of this paper is as follows. In section ^, we set up the general machinery 
we need, and prove the main theorems in the context of a discrete set 5 in a proper metric 
space. 

Section ^ explores this in the context of manifolds of constant curvature. The universal 
cover is M", S", or H", and the group G of deck transformations is a discrete group of isometrics 
(see [ doQ I). The discrete set S is the orbit of a point not fixed by any element of G under this 



discontinuous group. It is easy to see that the mediatrices in this case are totally geodesic 
spaces. From the basic property explained above, one can deduce that the n-th Brillouin zone 
is a fundamental domain for the group G va. X. 

In section ^, we calculate exactly the number of geodesies of length t that connect the 
origin to itself in two cases: the flat torus M'^/Z^ and the Riemann surfaces /T{p), for 
p G {2, 3, 5}. While these calculations could, of course, be done independent of our construction, 
the Brillouin zones help visualize the process. 

In the final section, we give a nontrivial example in the case of a non-Riemannian metric, 
and mention a connection to the question of how many integer solutions there are to the equation 

^fc _|_ ^A: _ £qj. £xed k. 

Acknowledgments: It is a pleasure to acknowledge useful conversations with Federico Bo- 
netto, Johann Dupont, Irwin Kra, Bernie Maskit, John Milnor, Chi-Han Sah, and Duncan 
Sands. Part of this work was carried out while Peter Veerman was visiting the Center for 
Physics and Biology at Rockefeller University and the Mathematics Department at SUNY 
Stony Brook; the authors are grateful for the hospitality of these institutions. 



2 Definitions and main results 



In this section, we prove that under very general conditions, Brillouin zones tile (as defined 
below) the space in which they are defined, generalizing an old result of Bieberbach [^]. With 
stronger assumptions, we prove that these tiles are in fact well-behaved sets: they are equal to 
the closure of their interior. 

Notation: Throughout this paper, we shall assume X is a path connected, proper (see below) 
metric space (with metric d{-, •)). We will make use the following notation: 

• Write an open r-neighborhood of a point xq as Nr{xo) = {x ^ X \ d{xo, x) < r}. 

• Define the circumference as Cr{xo) = {x £ X \ d{xo,x) = r}. 
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• Their union is the closed disk of radius r, denoted by Dr{xQ) = {x £ X \ d{xQ,x) < r}. 



Definition 2.1 A metric space X is proper if the distance function d{x, ■) is a proper map for 
every fixed x £ X. In particular, for every x £ X and r > 0, the closed ball Dr{x) is compact. 
Such a metric space is also sometimes called a geometry (See ^C(^). 

Note if X is proper, it is locally compact and complete. The converse also holds if X is a 



geodesic metric space (see Thm. 1.10 of |Gr|). The metric spaces considered here need not be 
geodesic. 

Definition 2.2 The space X is called metrically consistent if, for all x in X, R > r > 

in M, and for each a £ Cr{x), there is a z £ Cr{x) satisfying A'^^{^,a) (-2) ^ Xr{x) and 
Cd{z,a){z)r]CR{x) ={a}. 

This property is satisfied for any Riemannian metric. 

Note that any mediatrix La,fe separates X, that is: X — Lab contains at least two compo- 
nents (one containing the point a and the other b). 

Definition 2.3 We say that the mediatrix Lab is minimally separating if for any subset L CZ Lab 
with L 7^ Lab, set X — L has one component. 

If a mediatrix L it is minimally separating, then X — L has exactly two components. Note that 
if a separating set L <Z X contains a non-empty open set V, it cannot be minimally separating. 
For if A and B are disjoint open sets containing X — L, then so are A = A — {A Ci V) and 
B = B — (BnV). Now let x be any point in V. Then it is easy to see that the disjoint open sets 
AuV and B cover X — (L — x). Thus L — x separates X, so L could not have been minimal. 
We define the following sets 

= {x G X I d(0, x) - d{a, x) < 0} , L+ = {x G X | d(0, x) - d{a, x) > 0} . 



Definition 2.4 Two minimally separating sets L^a end L^b o-re topologically transversal if they 
are disjoint, or if for each x £ Loa H L^b and every neighborhood V of x, the sets Lq^ n L^^^ H V, 
Lq^ n n V, Lq^ n L^^^ n V, Lq^ n Lq^ n V are all nonempty. 

Usually, there will be a discrete set of points S = {xi}-^j in X which will be of interest. 
By discrete we mean that any compact subset of X contains finitely many points of S. Note 
that if liminf d{a, b) > 0, then S is discrete. 

a,beS 

Definition 2.5 We say a proper, path connected metric space X is Brillouin if it satisfies the 
following conditions: 

1: X is metrically consistent. 

2: For all a, b in X, the mediatrices Lab o'^e minimally separating sets. 
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3: For any three distinct points 0, a, and b in X, the mediatrices Loa o-nd Lq^ are topologically 
transversal. 

The last two conditions in the above definition may be weakened to apply only to those 
mediatrices Lab where a and b in S. In this case, we will say that X is Brillouin over S, if it is 
not obvious from the context. 




1 /| I ' Al 



Figure 2.1: The set L(o,o),{a,a) 
quarter-planes. 



contains two 





Figure 2.2: L(o,o),(4,6) (thin sohd Une) and Figure 2.3: The mediatrices Loa for with the Manhattan 
i(o,o),(2,4) (thick grey line) are not transverse, metric and a in the lattice { (m, nV2)} . 



Example 2.6 Equip with the "Manhattan metric", that is, d{p,q) = \pi — qi\ + \p2 — 92|- 
In this metric, a circle Cr{p) is a diamond of side length r\/2 centered at p, so condition 1 
is satisfied. However, condition 2 fails: if the coordinates of a point a are equal, then Loa 
consists of a line segment and two quarter-planes (see Fig. [2.11) . If the discrete set S contains 
no such points, we can still run into trouble with topological transversality. For example, the 
mediatrices -Z>{o,o),{2,4) and -Z^(o,o),(4,6) both contain the ray {(t, 1) | f > 4} (Fig. ^.21) . But, if we 
are careful, we can ensure that the space is Brillouin over S. To achieve this, if (0, 0) is the 
basepoint, we must have that for all pairs (ai,a2) and (6i,?*2) in ai — a2 ^ bi — 62- For 
example, take S to be an irrational lattice such as {(m, n\/2) | m,n G Z}. It is easy to check 
that for all a,b £ S, the properties of definition |2.5| are true. (From this example, we see that 
to do well in Manhattan, one should be carefully irrational.) 

As mentioned in the introduction, for each xq G 5*, the mediatrices L^^a give a partition 
of X. Informally, those elements of the partition which are reached by crossing n — 1 mediatrices 
from Xq form the n-th Brillouin zone, Bn{xo). This definition is impractical, in part because 
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a path may cross several mediatrices simultaneously, or the same mediatrix more than once. 
Instead, we will use a definition given in terms of the number of elements of S which are nearest 



to X. This definition is equivalent to the informal one when X is Brillouin over S (see Prop. 2.13 
below). We use the notation #{S) to denote the cardinality of the set S. 





Figure 2.4: Here we illustrate the definition of the sets bn{xo) and Bn{xo) for the lattice in R^. In both 
pictures, the circle Cd(x,xo)(x) is drawn, and the basepoint xq lies in the center of the square at the lower left. 
On the left side, the point x (marked by a small cross) lies in 65, and # {Nr{x) n S) = 4, while xq is the only 
point of S on the circle. On the right, we have m — 4 and ^ = 8, so a; lies in all of the sets B^,Be,, . . . , B\2- 



Definition 2.7 Let x £ X, let n be a positive integer, n < #{S), and let r = d{x,XQ). Then 
define the sets bn{xQ) and Bn^xo) as follows. 

• X £ bn{xo) <^=^ # {Nr{x) H S) = u — 1 and Cr{x) n S = {xq}. 

• X £ Bn{xo) <S=^> {Nr{x) n S) = m and # {Cr{x) n S) = £, wherel,meZ^ 
with m + 1 < n < m + i. 

Here the point xq is called the base point, and the set i?„(rco) is the n-th Brillouin zone 
with base point xq. Note that in the second part, if m = n — 1 and i = 1, then x £ 6„(xo). 
So bn{xo) ^ Bn{xo). Note also that the complement of bn{xo) in Bn{xo) consists of subsets of 
mediatrices (see Def. [l.ll) . Note also that 6n(xo) is open and that Bn{xo) is closed. Finally, 
observe that for fixed xq the sets bn(xo) are disjoint, but the sets Bn{xo) are not. In what 
follows it will be proved that {-Bn(a;o)}„>o cover the space X. Thus we can assign to each point 
X its Brillouin index as the the largest n for which x G Bn{xo). This definition was first given 



m 



The following lemma, which follows immediately from Def. 2.7, explains a basic feature 
of the zones, namely that they are concentric in in a weak sense. This property is also apparent 
from the figures. 
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Lemma 2.8 Any continuous path from xq to Bn{xQ) intersects -B„_i(xo). 

The Brillouin zones actually form a covering of X by non-overlapping closed sets in various 
ways. This is proved in parts. The next two results assert that the zones B cover X, but the 
zones b do not. The first of these is an immediate consequence of the definitions. The second 
is more surprising and leads to corollary 3.4. The fact that the Bi{xn) is the closure of bi{xn) 
(and thus that the interiors do not overlap) is proved in proposition |2.12| . We will use the word 
"tiling" for a covering by non-overlapping closed sets. 



Lemma 2.9 For fixed n the Brillouin zones tile X in the following sense: 
[jBi{xn)=X and 6i(a;„) n 6j(x„) = ifi^j- 





Figure 2.5: This example illustrates Lemma ^^9] and Thm. |2.1C| . Let S be the discrete set {(m,0)} U 
{(0, n)} ,m,n G Z in the Euclidean plane. On the left is the tiling given by _Bi(0, 0) and in the middle is the 

The 



2.10 



tiling by _Bi(2,0). In both cases, &2 is shaded. On the right is the tiling given by B2{xi) as in Thm. 
sets 62(0, 0), 62(1, 0), and 62(2, 0) ha. ve been shaded. Note that this S does not correspond to a group, nor does it 
satisfy the hypotheses of Prop. 2.11, because there are no isometries which permute S and do not fix the origin. 



Theorem 2.10 Let X he a proper, path connected metric space let S = {xjj-gj he a discrete 
set. Then, for fixed n < #{S), the sets {Bno{xi)}-^j tile X in the following sense: 

[J Bn{xi) = X and bn{xi) n bn{xj) = ifi^ j- 

i 

Proof: First, we show that for any fixed n > and each x £ X, there is an Xj S 5 with 
X £ Bn{xi). Re-index S so that if = {xi,X2,X3, . . .} and i < j, then d{x,Xi) < d{x,Xj). This 
can be done; since 5 is a discrete subset and closed balls Dc{xi) are compact, the subsets of S 
with d{x,Xi) < c are all finite. Let rj = d{x,Xi). We will show that x E Bn{xn)- 

Note that r„ > r„_i. Suppose first that r„ > r„_i, then Nr„{x) n S contains exactly 
n — 1 points, and x„ G Cr„{x) Pi S. Thus x £ Bn{xn)- Note that if r„+i > r„, then we would 
have x £ bn{xn) C Bn{Xn)- 

If, on the other hand, r„ = r„_i, then there is a A; > so that r„ = r„_i = . . . = r„_fc, 
and so # (iVr„(x) f^S) = n- k- l<n-l. But then # (Cr„(x) n S") > A; + 1, and hence 
X £ Bn{xn) as desired. 
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For the second part, we show that hn{xi)r\bn{xj) = 0. If not, then there is a point x in their 
intersection. If r, = rj, then Xi = Xj, because by the definition of bn{xk), {xk} = Cr^{x) n S. 
If not, then rj < rj. In this case, Xi S Dt.^{x) C Nr-{x) . Thus, since # {Nr^{x) Ci S) = n — 1, 
Nr- (x) must contain at least n points of S, a contradiction. □ 

The next result indicates how this notion of tiling is related to the notion of fundamental 
domain. 



Proposition 2.11 Let S be a discrete set in a metric space X as in theorem 2.1(\ . Suppose 



that for each Xi in S there is an isometry gi of X such that gi{xQ) = xi, gi permutes S and 
the only gi which leaves xq fixed is the identity. Then there is a set F (the fundamental set), 
satisfying: 

bn{xo)^ F ^ Bn{xo) with 
\Jgi{F) = X and g,{F) r\ gj{F) = % [i ^ j). 



Proof: Suppose that x € 6n(xo). From definition |2.7| and the fact that the gi are isometrics, 
we see that this is equivalent to gi{x) £ bn{xi). Thus gi{bn{xQ)) = bn{xi). Now apply theorem 



2.1C . A similar reasoning proves the statement for Bn{xo). □ 



Remark: The fundamental set F is not necessarily connected. Also, note that it follows from 



this proposition that Bi{xo) is scissors congruent to Bj{xo) (see [^ah[| for a discussion of scissors 
congruence). In particular, this implies immediately that the Bi all have the same area. Note 
that this result does not hold if S is not generated by a group of isometrics. For example, 
the "face-centered cubic" structures common to many crystals do not satisfy this, and the 



corresponding B^ have different volumes, although they still tile the space. See also figure 2.5 



The preceding results are valid without the transversality conditions on the metric, as the 
reader can check. We will now use the additional properties of the metric to derive a regularity 
result. 

Proposition 2.12 If X is Brillouin (over S), then Bn{xi) is the closure ofbn{xi). 

Proof: Let x £ Bn{xo)- If x G bn{xo), we are done. So suppose x bn{xo), and then 
by Def. (with r = d{x,xo)), Nr{x) contains m < n — 1 points of S and Cr{x) contains 
I > n — m with i > 2 points of S, including xq. Thus x lies in the intersection of £ — 1 > 1 
mediatrices {ixo,xJ-Ii • 

Suppose i = 2, so X lies on a single mediatrix Li. Then by condition 2, x is in the 
closure of and Lf. Thus if x' G Lj~ and close enough to x, then x' e 6^+1(3^0)- If x' G 
and sufficiently close to x, then x' £ 6^+2 (3^0) • From our choices, we see that m = n — 1 or 
m = n — 2. In either case, we are done. 

The proof for £ > 2 is slightly more involved (we have to use transversality), but the idea 
is the same. The point x lies in the intersection of at least two mediatrices. Essentially, we 
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need show that the mediatrices disentangle in every neighborhood of x. We wih prove that for 
all j G {1, ...,£}, a; G hm+j{xQ). Note that n — £<m<n — 1, and so showing this is sufficient. 

Again, since mediatrices are minimally separating, we have that for all i G {1,...,^ — 1} 
there exist open sets and in any neighborhood of x for which 

d(x,Xo) < d{x,Xi) for x £ L~ and d{x,Xo) > d{x,Xi) for x € Lf. 

Noting that Xi G Cr{x), then from condition 1 of Def. |2.7| , we have that for any small p > 0, 
there is a point z~ £ Cp{x) for which d{z^,XQ) < d{z~,Xi) for alH G {1, . . . , £}. Thus Hi 
is nonempty, and is contained in bm+iixo)- By the same principle, we can choose a point 
G Cr{x) for which d(z+,xo) > d{z^ ,Xi). As a consequence, C\iLf is nonempty, and is a 
subset of hm+i{xQ). 

Now, 6^+1(2^0) and &m+^(a;o) are disjoint. If £ > 1, their complement in a neighborhood of 
X must contain at least two mediatrices. By transversality (condition 3), these cannot coincide. 
Thus V must be intersected by some 6^_|_2(xo) for i G -[1, . . . By induction one finishes the 
argument. □ 



Remark: The references to mediatrices in the above proof may make it appear that we are 
relying on the informal definition of Bn rather than Definition 2/7. In fact, we are only using 
the fact that Bn \ bn consists of subsets of mediatrices, which was noted following Def. |2.7| . In 
fact, for Brillouin spaces, the informal definition of 6„ by mediatrices and Definition |2.7| are 
equivalent, as the following proposition shows. 



Proposition 2.13 Let X be Brillouin over S. Then there is a path 7 from x to xq which 
crosses exactly n — 1 mediatrices with no multiple crossings if and only if x £ 6„(a;o). By "no 
multiple crossings" we mean that if Li and Lj are distinct mediatrices, 7 intersects Li in at 
most a single point, and 'j D Li ^ 'j f] Lj. 



Proof: First, suppose there is such a path 7 from xq to x. As long as 7(t) crosses no mediatrices, 
the number of elements of S contained in Dr{'y{t)) remains constant. If this number is i, then 
7(i) G 5i(xo). Each time j{t) crosses a single mediatrix, this number increases by 1. Thus, 
after crossing n — 1 mediatrices, 7(i) will be in 6„(xo). 

Now, let X G bnixo). We must construct a path from x to xq which crosses exactly n — 1 
mediatrices. Let 6* denote the connected component of bn{xo) which contains x, and choose a 
point y in i?„_i(xo) H 96*. Such a point always exists by Lemma 2.8 and Prop. 2.12. In fact, 
because dbn consists of subsets of mediatrices, topological transversality of mediatrices allows 
us to choose y so that it lies on a single mediatrix. 

There certainly is a path connecting x to y which crosses no mediatrices; by similar 
reasoning, we can choose x' G bn-i{x()) and a path from x' to y which crosses no mediatrices. 
We have now constructed a path from x £ bnto x' € which crosses exactly one mediatrix — 
repeating the argument n — 2 times finishes the proof. □ 
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3 Brillouin zones in spaces of constant curvature 

In this section X will be assumed to be one of M", S", or H", all equipped with the standard 
metric, and let G be a discontinuous group of isometrics of X. Denote the quotient X/G with 
the induced metric by (M, g). Then the construction of lifting to the universal cover, as outlined 
in the introduction, applies naturally to (M, g). In this section we describe focusing of geodesies 
in {M,g) by Brillouin zones in X. The discrete set S is given by the orbit of a chosen point in 
X (which we will call the origin) under the group of deck-transformations G. The fact that the 



Brillouin zones are fundamental domains is now a direct corollary of proposition 2.11 



The regularity conditions of Def. 2^ are easily verified in the present context. We do this 

first. 

Lemma 3.1 If X is either M."-', §"", orH", then a mediatrix Lab in X is an (n—l)- dimensional, 
totally geodesic subspace consisting of one component, and X — L^t has two components. 

Proof: This is easy to see if we change coordinates by an isometry of X, putting a and b 
in a convenient position, say as x and —x. The mediatrix L^-x is easily seen to satisfy the 
conditions (in the case of it is the equator, and for the others, it is a hyperplane). The 
conclusion follows. □ 



Proposition 3.2 All such spaces X are Brillouin (see definition 
Proof: As remarked before, the first condition is satisfied for any Riemannian metric. The 



are 



second condition is also easy. It suffices to observe that the subspaces of Lemma |3.1 
minimally separating. 

To prove the topological transversality condition, note that lemma 3.1 implies that if Loa 



and Lob coincide in some open set, then they are equal. One easily sees that then we must have 
Loa = Lqij = Lab = L. Applying the second part of the lemma to Loa, Lqi,-, and Lab-, we see that 
they must separate X in at least 3 components, namely one for each point a, 6, and 0. This 
leads to a contradiction. □ 

Remark: Note that in fact, the mediatrices are transversal in the usual sense. If Loa and L^b 
coincide in an open set, then their tangent spaces also coincide at some point. Uniqueness of 
solutions of second order differential equations then implies Loa = L^b- 

Recall that a metric space X is called rigid if the only isometry which fixes each point 
of a nonempty open subset of X is the identity. It is not hard to see that S" , H", and M" are 



rigid spaces. See [Ra] for more details of rigid metric spaces and for the proof of the following 
result. Recall that the stabilizer in G of a point x ^ X consists of those elements of G that fix 

X. 

Proposition 3.3 Let G be a discontinuous group of isometrics of a rigid metric space X. Then 
there exists a point y of X whose stabilizer Gy is consists of the identity. 
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We now return to Brillouin zones as defined in the last section. Recall that G is a group 
of isometries of X that acts discontinuously on points in X. Let xq be a point in X whose 
stabilizer under the action of G is trivial. Let x G X be a point in the orbit of xq under G, that 
is to say, x = 7(xo) for some 7 G G — {e}, where e is the identity of G. Let [xq, x\ be a geodesic 
segment of minimal length whose endpoints are xq and x. Then Bn^xo), the n-th Brillouin 
zone relative to xq , is the set of points in X such that the segment [xq , x] intercepts exactly 
n — 1 mediatrices Lx^^y, where y is in the orbit of xq under the group G. Proposition ^.11 
immediately implies the most important fact about Brillouin zones in this setting. 





Figure 3.1: Brillouin zones for PSL{2,'L) in the hyperbolic disk. We have transported the "usual" upper 
half- plane representation using the map z 1— > On the left is shown the sets Bn{j), which give fundamental 

On the right, is taken as a basepoint. Since the origin has a non-trivial stabilizer, the 



domains as in Cor. 3.4 



corresponding Brillouin zones give a double cover of the fundamental domains. 



Corollary 3.4 Let X he R", S", or H", and let G he a discontinuous group of isometries of 
X . Let Xq £ X he such that its stabilizer Gxq under G is trivial. Then for every positive integer 
n, the n-th Brillouin zone i?„(xo) is a fundamental domain for the action of G on points in X . 
Its boundary is the union of pieces of totally geodesic subspaces and equals the boundary of its 
interior. 

Remark: The first Brillouin zone -Bi(O) is the usual Dirichlet fundamental domain for the action 
of G. Furthermore, even when Gxq is not trivial, B^ixo) is a A:-fold cover of a fundamental 
domain. 

As pointed out in the introduction, the number of geodesies that focus in a certain point 
is counted in the lift. So if a given point x £ X is intersected by n mediatrices, it is reached by 
n + 1 geodesies of length d{0,x) emanating from the reference point (the origin). In the next 
section, we give more specific examples of this. 

Finally, we state a conjecture. 
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Conjecture 3.5 Let {X,g) be the universal cover of a d- dimensional smooth Riemannian man- 
ifold (M, g) as described in the construction. For a generic metric g on M , no more than d 
mediatrices intersect in any given point y of X. 

This conjecture acquires perhaps even more interest (and certainly more structure), when 
one restricts the collection of metrics on M to conformal ones ( [ [Mas| ] ) . A result in this direction 
for M = R^/Z^ can be found in poll] . 



4 Focusing in two Riemannian examples 

In this section, we give two examples (one of them new as far as we know) of focusing. Suppose 
that at t = geodesies start emanating in all possible directions from the a point. At certain 
times ti,t2,----, we may see geodesies returning to that point. We derive expressions for the 
number of geodesies returning at t„ in two cases. First, as introductory example we will discuss 
this for M = R^/Z^ (a more complete discussion of this example can be found in jPisj ). 
Second, we will deal with a much more unusual example, namely M = IHI^/r(A;), where r(A:) is 
a subgroup of PSL{2,Z) called the principal congruence subgroup of level k (defined in more 
detail below). We note that it seems to be considerably harder to count geodesies that focus in 
points other than our basepoint. 

Before continuing, consider the classical problem of counting Rg{n), the number of solu- 
tions in of 

2 , 2 

p + q = n. 

Let 

k £ 

i=i j=i 

be the prime decomposition of the number n, where pi = l(mod4) and qi = 3(mod4). The 
following classical result of Gauss (see |NZ^) will be very useful. 



Lemma 4.1 Rg{n) is zero whenever n is not an integer, or any of the ji is odd. Otherwise, 

k 

1=1 

Example 4.2 Choose an origin in M = R^/Z^ and lift it to the origin in R . Our discrete set 
S is then Z^. Let Px{t) be the number of geodesies of length t that connect the origin to the 
point X S M. 

Proposition 4.3 In the flat torus, the number of geodesies of length t that connect the the 
origin to itself is given by 



Po{t) 



Rg{t^), ift^en 

0, otherwise. 
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Proof: Notice that by definition geodesies of length t leaving from the origin in reach the 
points contained in Ct{0). Only if t"^ is an integer does this circle intersect points in I?. □ 



Example 4.4 We now turn to the next example. Recall that PSL{2, Z) can be identified 
with the group of two by two matrices with integer entries and determinant one, and with 
multiplication by —1 as equivalence. For each k, the group T{k) is the subgroup of PSL{2,'L) 
given by 

a = d= 1 (mod A;), 6 = c = (mod k)^ . 

This group has important applications in number theory. The action of r(A;) on is given by 
the Mobius transformations 

, . az + b , / a b \ . 

q(z] = where , G 1 (fc . 

^ ' cz + d \ c d J ^ ' 

We point out that for k = 2, 3, or 5, the surface /T{k) is a sphere with 3, 4, or 12 punctures 
(see il]). 

We will find it more convenient to work in the hyperbolic disk D'^, which is the universal 
cover of EI^/r(A;). We shall choose a representation of r(/c) in the disk so that i E corresponds 
to the origin. This will allow us to determine the focusing of the geodesies which emanate from 
i. Note that the surface EI^/r(fc) has special symmetries with respect to i: for example, i is the 
unique point fixed by the order 2 element of PSL[2,'L). 




m 



G PSL{2,Z) 



Figure 4.1: The orbit of i under r(2) transported to the hyperbohc disk, and the corresponding Brillouin 
zones. Each zone Bn forms a fundamental domain for a 3-punctured sphere. 
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Lemma 4.5 The action of the fundamental group of the surface /T{k) can he represented 



as 



r — IS p + iq 
p — iq r + is 



r+p=l (mod k), r -p=l (mod k) ,2,2,1 2,2 
, -n) A A -n) A i,\ andp^ + q^ + 1 = r^ + s^ 

s + q = U (mod fej, s — q = U (mod fcj 



acting on B^. We shall denote this particular representation as the group rQ{k). 



Proof: Following the conventions in |Be|, define 



(h-.B^^M^, (hiz) = i . 

, ^\ ; -z + 1 

Push back the transformation g G ^{k) from to by g' — > (j)~'^g(f) to obtain a representation 
of 5 G T{k) as a transformation acting on D^. The matrix representation of this transformation 
is given by: 

a+d I ■ b—c a—d ■ b+c 

2 ~^ 2 2 2 

An 



a—d I ■ b+c a+c 
2 ' 2 2 



■.b-c 



where det A„ = 1, since this matrix is conjugate to g, whose determinant is equal to 1. Let 



p={a- d)/2 
r = {a + d)/2 



q=-{b + c)/2 
s = -(h-c)/2 



and Ag now written as 



A, 



r — is p + iq 
p — iq r + is 



Here the numbers p, q, r, s are in Z and must satisfy the following congruence conditions: 

r + p = 1 (mod k), r — p = 1 (mod k) 
s + q = (mod k), s — q = (mod k) 

Since the determinant of Ag is equal to 1, we must also have 

p +g + l = r +s. 



□ 



We need another auxiliary result before we state the main result of this section. 

Lemma 4.6 Let {p, q) and (r, s) be two points in 1? such that the integers A = p^ + q^ and 
B = r"^ + s"^ are relatively prime, and let be a rotation. Now ^p{p,q) = {p',q') and ip{r,s) = 
{r',s') are in 1? if and only if (p is a rotation by a multiple of 71/2. 
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Proof: Let c be the cosine of the angle of rotation. We have 

p'p + q'q r'r + s's 



A 



B 



Thus if p'p + q'q and r'r + s's are not both equal to zero, 

p'p + q'q A 
r'r + s's B 

Because A and B are relatively prime and surely \p'p + q'q\ is less than or equal to A, and 
similarly for B, we have that 



p'p + q'q = ±A and r'r + s's = ±B. 



This implies the result. 



□ 



Now we define a counter just as before. Choose a lift of M = B^/r0(A;) so that G M 
lifts to G . Let (*) be the number of geodesies of length t that connect the origin to the 
point X G M. 



512 
256 
128 
64 
32 
16 
8 
4 



—I 1 1 1 1 1 1 1 1 1 r- 



—I 1 1 1 1 1— 



_i I I I I I I I I I i_ 



_i I I I I i_ 



512 ■ 
256 ■ 
128 ■ 
64 ■ 
32 ■ 
16 ■ 
8 ■ 
4 ■ 



1 1 1 1 1 1 1 1— 



J I I I I I I l_ 



012345678 9 10 11 12 13 14 15 16 17 18 19 2021 22232425 



012345678 



Figure 4.2: The non-zero values of po{t) for t < 25 (left) and 7o(t) for t < 8 (right), which count how many 
geodesies of length t connect the origin to itself in the R^/Z^ and n'^/rQ{2), respectively. 



Theorem 4.7 In the surface E[^/r(/c), the number of geodesies which connect the point i G 
to itself is given by 

iRg (cosh^ t - l)Rg (cosh^ t) fork = 2 
\Rg (^coshjt^^ (cosh2 1) fork = 3 

iRg ^p(cosh2 1) fork = 5 

Note that in all cases, the number is nonzero only if cosh^ t G N. 
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Proof: We shall work in the disk, rather than in H^. Let S be the orbit of G under r0(A:). 
Then the number of such geodesies is exactly the number of distinct points of S which lie on 
on the circle Ct{0) of radius t and centered at the origin. 

If X G S, then by Lemma 4.5, it is of the form with p, q, r, s integers satisfying 
+ = + — 1. Let n be their common value, that is, 

n = p'^ + q'^ = r'^ + s"^ - 1. 

We will first count the number of 4-tuples {p, q, r, s) that solve this equation, momentarily 
ignoring the congruence conditions. 

Note that the point x has Euclidean distance to the origin given by 

. ,2 p^ + q^ n 
\x\ — — 



+ n + 1 

The hyperbolic length of the geodesic connecting the origin to this point is 

t = arctanh (|a;|e). 



This gives that 7o(i) is only non-zero when t = arctanh y o'^! equivalently, n = cosh^ t — 1. 

To count the number of intersections of Ct(0) with S for these values of t, observe that 
we can use Gauss' result to count the number of pairs {p, q) such that p^ + q^ = n. This number 
is given by Rg{n). For each such pair (p, g), we have a number of choices to form 

p + iq 

X = 

r + is 

By the above, this number is equal to Rg{n + 1). Thus, 70 (t) is at most Rg{n)Rg{n + 1). 

However, we have over-counted: some of our choices for p, g, r, s represent the same point 
2; E S*, and some of them may not satisfy the congruence conditions, which we have so far 
ignored. We will first account for the multiple representations, and then account for the con- 
gruence relations. 

Let p,q,r,s G Z be as above, giving a point x = at distance t = arctanh \J 

from the origin. If we multiply the numerator and denominator of x by e*^, then x will remain 
unchanged. Because of the requirement that + = + — 1 = n, this is the only invariant, 
and by Lemma 4.6, 9 must be a multiple of ^ for the numerator and denominator to remain 



Gaussian integers. We see that in our counting, we have represented our point x in 4 different 
ways: 

_ p + iq _ —q + ip_ —p — iq _ q — ip 
r + is —s+ir —r — is s — ir 

meaning we have over-counted by a factor of at least 4. 
Now we account for the congruence conditions. 

First, consider the case k = 2. Note that q + s = (mod 2) if and only if p + r = 
1 (mod 2), because p^ + + 1 = + s^, so we need only check this one condition. If the 
representation fails to satisfy our parity condition, then q+s = 1 ( mod 2) and consequently 
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p + r = ( mod 2). This means that the representation of this same point does satisfy the 

parity conditions, giving exactly \Rg{cosb? t — l)i?g(cosh^ t) distinct points of S at distance t 
from the origin. 

If A; = 3, then since k is odd, the congruence conditions on p, g, r, and s imply that 
r = l(mod 3) and p = q = s = O(mod 3). 

Note that the equation 

p^ + = n and p = g = (mod 3) 

wiU be satisfied exactly Rg{n/3'^) times. (Recall that if n/9 is not an integer, Rg{n/9) = 0.) 

For fixed n, let (p, q) be any one of the solutions. We need to decide how many solutions 
the equation 



r 



^ + = n + 1 with r = 1 (mod 3) and s = (mod 3) 



admits. The solution of the first equation implies that 3 divides n. Thus + = 1 (mod 3). 
Consequently, we have 4 choices mod3 for the pair (r, s), namely (0, 1), (1, 0), (0, 2), and (2, 0). 

Let {p, q, r, s) G be any solution to n = p^+q^ = r^+s^— 1 withp = g = 0( mod 3). 

For each choice of {p, q), we have exactly Rg{n+1) choices of (r, s). Now let R denote the product 
of the rotations by 7r/2 on each of the components of x Z^. Using Lemma 4.6, we see that all 
such solutions can be obtained from just one by applying R repeatedly. It is easy to check that 
each quadruple of solutions thus constructed runs exactly once through the above list. Since 
precisely one out of the four associated solutions is compatible with the conditions, the total 
number of solutions is exactly: 

Using the relationship between the Euclidean distance and the Poincare length as before gives 
the result. 

If A; = 5, the proof for A; = 3 can be literally transcribed to obtain the result. □ 

Remark: Note that the above results do not hold if k is not one of the cases mentioned. The 
primary difficulty is that for prime k >7, there are solutions which are not related by applying 
the rotation R. However, the argument does give iRg ((cosh^ t - 1)/A;2) Rg{cosh'^ t) as an upper 
bound for BI^/r(A:) when k is an odd prime. Note that the surface E[^/r(A;) is of genus if and 
only if A: < 5 (see iFI 



5 Non-Riemannian examples 



The present context is certainly not restricted to Riemannian metrics. As an indicator of this 
we now discuss a different set of examples. 
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Let k he a positive number greater than one. Equip with the distance function 



X -y 



xi - yi\^ + \x2 - y2f) 



l/k 



and let the discrete set S be given by Z^. For k not equal to 2, this is not a Riemannian metric, 
yet all conclusions of section |2| hold. In particular, each Brillouin zone forms a fundamental 
domain. Note that determining the zones by inspecting the picture requires close attention! 





Figure 5.1: Brillouin zones for the lattice 7? in with the metric (|x-i — yif + \x2 — J/2 
Figure 2.3 and Example 2.6, which deal with the case fe = 1, the "Manhattan metric". 



4\l/4 



See also 



Now the problem of determining Ct{0) Ci S for any given t is unsolved for general k. In 
fact, even for certain integer values of k greater than 2, it is not known whether Ct(0) n S ever 
contains at least two points that are not related by the symmetries of the problem. For /c = 4, 
the smallest t for which Ct{0) n S has at least two (unrelated) solutions is given by 



133'^ + 134"^ 



158^ + 59^ 



However, for fc > 5, it unknown whether this can happen at all (see [ pW| ]). 

There are some things that can be said, however. In the situation where k £ {3, 4, 5, . . .}, 
the mediatrices intersect the coordinate axes only in irrational points or in multiples of 1/2. 
For if X = {p/q, 0) is a point of a mediatrix -^^(o.o), (01,02)) '^^ have 

\p\^ = \p - qai\^ + \a2q\^ (p/0,g/0) . 

By Fermat's Last Theorem, this has no solution unless either p = qai or 02 = 0. In the first 
case, I = ±02, which can only occur if the lattice point is of the form (02, ±02)- If ^2 = 0, then 
- = In particular, there is no nontrivial focusing along the axes. 
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To compute Figure 5.1, we took advantage of the smoothness of the metric. Not all 
metrics are sufficiently smooth for this procedure to work. Even for Riemannian metrics, in 
general the distance function is only Lipschitz, which will not be sufficiently smooth. 

For each a = (01,02) G Z^, define a Hamiltonian: 

Ha{x) =11 X — a II — II x II . 

The mediatrix Loa corresponds to the level set Ha{x) = 0. Because Ha{x) is smooth, we have 
uniqueness of solutions to Hamilton's equations. In the current situation, where the dimension 
is two, the level set consists of one orbit. Thus, one can produce the mediatrix by numerically 
tracing the zero energy orbits of the above Hamiltonian. 

As mentioned above, for a general Riemannian metric, the distance function is only 
Lipschitz. This means we have no guarantee that the solutions of the above differential equation 
are unique. Indeed, there are examples of multiply connected Riemannian manifolds with self- 
intersecting mediatrices, as will be shown in a forthcoming work. 
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